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ABSTRACT 



' Context. The characterization of circumstellar disks around young stellar objects can provide important information about the process 

|_ of star formation and the possible formation of planetary systems. 

Aims. We investigate the spatial structure and the spectral energy distribution of a newly discovered edge-on circumstellar disk around 
an optically invisible young stellar object that is embedded in a dark cloud in the Carina Nebula. 

Methods. The disk object was serendipitously discovered in our deep near-IR imaging survey of the Carina Nebula obtained with 



o 



HAWK-I at the ESO VLT. Whereas the object was detected as an apparently point-like source in earlier infrared observations, only 
the superb image quality (FWHM « 0.5") of the HAWK-I data could reveal, for the first time, the peculiar morphology of the object. 
It consists of a very red point-like central source that is surrounded by a roughly spherical nebula, which is intersected by a remarkable 
dark lane through the center. We construct the spectral energy distribution of the object from 1 /jm to 870 fjm and perform a detailed 
[ radiative transfer modeling of the spectral energy distribution and the source morphology. 

\ > j" > ) . Results. The observed object morphology in the near-IR images clearly suggests a young stellar object that is embedded in an 

• extended, roughly spherical envelope and surrounded by a large circumstellar disk with a diameter of as 5500 AU that is seen nearly 

£ — ' edge-on. The radiative transfer modeling shows that the central object is highly luminous and thus must be a massive young stellar 

£>f) \ object, most likely in the range 10 - 15 M B . The circumstellar disk has a mass of about 2 M e . 

Conclusions. The disk object in Carina is one of the most massive young stellar objects for which a circumstellar disk has been 
detected so far. The size and mass of the disk are very large compared to the corresponding values found for most other similar 
objects. These results support the assumption that 10 - 15 M stars can form via accretion from a circumstellar disk. 

Key words. Stars: formation - Stars: circumstellar matter - Stars: pre-main sequence - ISM: individual objects: NGC 3372 
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1. Introduction planets may be forming, evolves with time and is finally dis- 
^ ■ persed (see Dullemond et al. 2007). 
During the past dozen years, a substantial number of circum- 
stellar accretion disks around young stellar objects have been So far ' nearl y a11 wel1 characterized circumstellar disks sur- 
spatially resolved with a number of different techniques such round low " or intermediate mass stellar objects, with less than 
as direct optical/infrared imaging, (sub-)mm or radio-mapping, about 5 - 8 M . Only a very limited number of disk detec- 
radio-interferometry, and (more recently) infrared interferome- tlons have been reported around objects with higher masses in 
try (see, e.g., Stapelfeldt et al. 1998; Zinneckeret al. 1999; Wolf the ran 8 e M ~ 8 - 20 M Q (see, e.g., Chini et al. 2004; Jiang 
et al. 2003; Leinert et al. 2004; Preibisch et al. 2006; Kraus et et al - 2005 ; Nielbock et al. 2007; Kraus et al. 2010; Davies et 
al. 2008; Ratzka et al. 2009; Gramajo et al. 2010; Quanz et al. al - 2010 )- However, in many of these cases the mass estimates 
2010). Spectroscopic observations have also provided important for the central star quite uncertain; the true mass may in 
kinematical evidence of accretion disks (see, e.g., Chandler et fact be considerably smaller than initially estimated (see, e.g., 
al. 1995; Bik & Thi 2004; Wheelwright et al. 2010; Davies et al. Sako et aL 2005 )- Another problem is that the observed flattened 
2010). These studies have greatly advanced our understanding structure around these more massive young stellar objects are 
of the physics of young stellar objects (YSOs), how the young not necessarily accretion disks, from which material is directly 
stars gain their final mass, how the circumstellar matter, in which transferred to the central star; they may rather represent a thick 

torus that could be the remnant of a flattened envelope (Allen et 
al. 2003). Considering even higher stellar masses, M ^ 20 M Q , 



* The near-infrared observations for this project were collected with no clear observational evidence for accretion disks has yet been 
the HAWK-I instrument on the VLT at Paranal Observatory, Chile, un- found (see, e.g., Cesaroni et al. 2007). This rarity and lack of 
der ESO program 60.A-9284(K). convincing disk detections is an important aspect in the long- 
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standing discussion about how massive stars form. A fundamen- 
tal problem is that the strong radiation pressure resulting from 
the very high luminosity of a massive protostars tends to halt 
the accretion flow and should thus limit the final stellar mass, at 
least in the case of spherical accretion (Kahn 1974). More recent 
theoretical results alleviate this problem by mechanisms such 
as non-spherical accretion, the "flash-light effect", and "photon- 
bubbles" (see, e.g., Yorke & Sonnhalter 2002; Krumholz et al. 
2009). The more accurate treatment of frequency-dependent ra- 
diative feedback by Kuiper et al. (2010) suggests that stars with 
masses well beyond the upper mass limit of spherical accretion 
can be formed by accretion. However, due to the limited numer- 
ical resolution of these simulations (e.g., 1 .27 AU in the study 
of Kuiper et al. 2010), the reliability of the numerical results is 
not fully clear. As discussed in Zinnecker & Yorke (2007), there 
are still strong indications that the formation of massive stars is 
not simply a scaled-up version of the low-mass star formation 
process. Alternative models for the formation of massive stars 
highlight the importance of the cluster environment during early 
star formation stages (see, e.g., Bonnell et al. 2004; Krumholz et 
al. 2010), and processes such as protostellar collisions and merg- 
ers may be possible in the central regions of very dense clusters 
(see, e.g., Bonnell et al. 1998; Bonnell & Bate 2005; Bally & 
Zinnecker 2005; Davies et al. 2006). 

In the context of these open questions, every new detec- 
tion and good characterization of a circumstellar disk around a 
massive (M S 8 M Q ) young stellar object provides an impor- 
tant new mosaic stone that helps to solve the puzzle of massive 
star formation. However, observations of massive young stellar 
objects are strongly hampered by the relative rarity of massive 
stars, their correspondingly larger average distances (as com- 
pared to low-mass stars), and the extremely short timescales on 
which massive protostars evolve. The expected lifetime of a disk 
around a massive young stellar object is very small, ^ 10 5 yrs, 
since the circumstellar material is quickly dispersed by the enor- 
mous luminosity and wind power of the young high-mass stars. 
Furthermore, due to their extreme youth, massive YSOs with 
disks are still deeply embedded in their natal cloud and thus op- 
tically invisible. 

The Great Nebula in Carina (NGC 3372; see Smith & Brooks 
2008, for an overview) is a superb location in which to study 
the physics of massive star formation. At a distance of 2.3 kpc, 
it represents the nearest southern region with a large popula- 
tion of massive stars, among them several of the most massive 
(M ~ 100 M ) and luminous stars known in our Galaxy. Most of 
these very massive stars reside in the clusters Trumpler 14 and 
16, which are thought to have ages of a few Myr (Tapia et al. 
2003). Recent sensitive infrared, sub-mm, and radio observa- 
tions showed clearly that the Carina Nebula Complex (CNC) is 
a site of ongoing star formation. The region contains more than 
10 5 M o of gas and dust (see Smith & Brooks 2008; Preibisch 
et al. 2011a). Several very young stellar objects (Megeath et al. 
1996; Mottram et al. 2007) and a spectacular young embedded 
cluster (the "Treasure Chest Cluster"; see Smith et al. 2005) have 
been found in the molecular clouds, a deep HST Ha imaging sur- 
vey revealed dozens of jet-driving young stellar objects (Smith 
et al. 2010a), and Spitzer surveys located numerous embedded 
protostars throughout the Carina complex (Smith et al. 2010b; 
Povich et al. 2011). A recent deep wide-field X-ray survey re- 
vealed at least ~ 1 1 000 young stars in this area (Townsley et al. 
201 1; Preibisch et al. 201 lb). The CNC is thus an ideal target in 
which to search for very young massive stars. 

We have recently used the ESO 8 m Very Large Telescope to 
perform the largest and deepest near-infrared survey of the CNC 



obtained so far. In our images, we discovered an optically invis- 
ible infrared source with a very interesting morphology, which 
seems to be a young stellar object surrounded by a huge cir- 
cumstellar disk seen nearly edge-on. It will be referred to as the 
"disk object" in the following text. In this paper we first discuss 
the morphology of the object and its environment, and then de- 
termine its spectral energy distribution from the near-infrared to 
the sub-mm range. We compare the spectral energy distribution 
to models of young stellar objects surrounded by circumstellar 
disks and envelopes. Finally, we present the results of our com- 
prehensive and detailed radiative transfer simulations in which 
we simultaneously modeled the spectral energy distribution and 
the observed near-infrared morphology of the object in order to 
determine basic parameters of the central young stellar object 
and the physical properties of the surrounding disk and enve- 
lope. 

2. Observations 

The near-infrared data presented in this paper were obtained 
in the context of our large survey of the Carina Nebula (see 
Preibisch et al. 2011c, for more details) with the instrument 
HAWK-I (see Kissler-Patig et al. 2008) at the ESO 8m Very 
Large Telescope. HAWK-I is equipped with a mosaic of four 
Hawaii 2RG 2048 x 2048 pixel detectors with a scale of 
0.106" per pixel. All data were processed and calibrated by the 
Cambridge Astronomical Survey Unit using pipeline software 
originally designed for the analysis of the UKIRT Infrared Deep 
Sky Survey. In the final images, objects as faint as / ~ 23 mag, 
H ~ 22 mag, and K s ~ 21 mag are detected at S /N = 3. 

The parts of the HAWK-I mosaic in which the disk object 
is located were obtained in the night of 27 January 2008. The 
image quality of these specific frames is particularly good; we 
measured the PSF size for several point-like sources near the 
disk object and found FWHM values of 0.37" in the /T s -band, 
0.45" in the //-band, and 0.52" in the /-band. This excellent im- 
age quality (for ground-based seeing-limited observations) was 
essential for the detection and morphological characterization of 
the disk object. 

In order to get more comprehensive information about the 
disk object we also analyzed archival data from the Spitzer Space 
Telescope (IRAC and MIPS maps) and optical Hubble Space 
Telescope (HST) images taken with WFPC2 and ACS. 

Finally, we also used our own deep sub-mm map of the CNC 
(Preibisch et al. 201 la) that we recently obtained with the Large 
APEX Bolometer Camera LABOCA at the APEX telescope. 
This map traces the 870//m emission at 18" angular resolution 
(= 0.2 pc at the distance of the Carina Nebula). 

3. Discovery of the disk object in the dark cloud 
near Cr 232 

The disk object was discovered in our visual inspection of the 
HAWK-I images due to its peculiar morphology. It is the bright- 
est member of a group of very red, deeply embedded objects that 
are surrounded by extended nebulosity, as shown in Fig. 1. The 
disk object is located at J2000 coordinates a = 10 h 44 m 31.06 s , 
5 = -59° 33' 10.5" and is seen as an extended, roughly spherical 
nebulosity with a diameter of about 4". The nebulosity is inter- 
sected by a prominent dark lane that divides it into two similarly 
large hemispheres. The dark lane can be clearly traced over a dis- 
tance of 3", corresponding to about 7000 AU. The global mor- 
phology of this object is similar to the famous 'Butterfly Nebula' 
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Fig. 1. Upper image: Optical image of the central part of the 
Carina Nebula, reproduced from the ESO photo release 1031 
(http://www.eso.org/public/images/etamosaicnm2/). North is up and 
east to the left. The box shows the area covered by the near-infrared 
HAWK-I image below. 

Lower image: RGB composite image constructed from the J- (blue), 
H- (green), and A",-band (red) HAWK-I images of the area around the 
disk object. (Note: the horizontal and vertical dark streaks are artifacts 
related to the dither pattern and the mosaicing process.) 



in Taurus, (a circumstellar disk surrounding a T Tauri star, see 
Wolf et al. 2003) or the 'Silhouette Disk' around a young star 
in M 17 (Chini et al. 2004), and is discussed in more detail in 
Sect. 4. 

3. 1 . The environment of the disk object 

In the optical images of this region (see Fig. 1) the disk object is 
completely invisible since it is embedded in a dark cloud in the 
region of the clustering Collinder 232. In order to show the struc- 
ture of this cloud, we plot in Fig. 2 the contours of the 870 /im 
emission (tracing the cold dust in the cloud) on the 8 /im Spitzer 
image (tracing the surface of the cloud). The south-eastern edge 
of the cloud is very sharp, straight and bright; its orientation is 
nearly perpendicular to the direction towards i] Car. This sug- 




Fig. 2. Negative grayscale representation of the 8 //m Spitzer image of 
the cloud in which the disk object is embedded, with overplotted con- 
tours of the 870 nm LABOCA map. The field-of-view is « 4.5' x 3.8'; 
north is up and east to the left. The sub-mm contour levels increase 
from 0.06 Jy/beam to 1.5 Jy/beam in equal steps. The position of the 
disk object is marked by the white circle. 



gests that this cloud edge is an ionization front caused by the 
strong UV irradiation from 77 Car and the O-type stars in Tr 16. 

The sub-mm emission of this cloud shows two peaks. The 
southern peak agrees very well with the position of the disk ob- 
ject, suggesting that this object dominates the sub-mm emission 
of this peak. The second, northern, peak is found 19" east and 
23" north of the disk object. It coincides with the center of the 
optical dark cloud and a corresponding minimum in the 8 /urn 
surface brightness (see Fig. 2). No embedded infrared sources 
are seen near the center of this second sub-mm peak. 

3.2. Previous detections of the disk object as an unresolved 
point-like source 

Since the disk object is a relatively bright infrared source, it 
was detected in many previous infrared observations of this re- 
gion, but the details of its spatial structure always remained com- 
pletely unresolved. Only due to the superior image quality of our 
HAWK-I data that provide an angular resolution of ~ 0.5", its 
peculiar morphology could be resolved. 

Figure 3 shows a small part of the H-band HAWK-I image 
around the disk source in which the locations of other sources are 
marked. The HAWK-I images show a bright point-like source 
only 2.15" to the south-east of the disk object. Its PSF shows 
no indications of being extended and it is clearly separated from 
the nebulosity surrounding the disk object. We thus believe that 
this point source is an unrelated star. Based on its position of or = 
10 h 44 m 31.3 s 6 = -59° 31' 11" we will refer to it as star 31.3-11 
in the following. This star is optically visible and much brighter 
than the disk object in the /-band, but less bright than the disk 
object in the ^-band. These rather blue colors may suggest it to 
be a foreground star. 

Inspection of the 2MASS images shows that the disk object 
and this nearby point source appear there as an unresolved blend, 
listed as 2MASS-J10443122-59331 13 in the point source cata- 
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Fig. 3. H-band HAWK-I image of the area around the disk object. The 
small (green) circles with 2" diameter mark the positions of the objects 
listed in the 2MASS point source catalog. The 5" diameter (red) circles 
mark the locations of the two 10.4 /jm sources detected by Mottram et 
al. (2007), and the 18" diameter (red) circle marks the peak of the sub- 
mm emission as seen in our LAB OCA map. A grid of J2000 coordinates 
is shown. 
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Fig. 4. RGB composite image constructed from the K s -, H-, and /-band 
HAWK-I images. 



log. Its position is just between the disk object and the south- 
eastern point source (as expected). 

The disk object is associated with the MSX source 
G287.4779-00.5463, which appears clearly extended (FWHM 
~ 40") in the MSX images. The object was also detected as 
a faint diffuse nebulosity in the thermal-infrared study of the 
Carina nebula by Rathborne et al. (2002), where it is denoted as 
source "N4". These authors performed a careful analysis of the 
MSX and IRAS data with the aim to separate the flux from the 
compact infrared source from the surrounding diffuse emission 
and derived the spectral energy distribution of the source from 
8/um to 100/um (their Fig. 12). They modeled the object with 
a two-component blackbody fit to the SED and derived a total 
luminosity of 41 000 L e and a spectral type of 09-09.5 for the 
central object. This result was the first hint that the disk object is 
actually a massive young stellar object. 

Mottram et al. (2007) performed 10.4 /im mid-infrared ob- 
servations of G287.4779-00.5463. Their images were obtained 
with TIMMI2 at the 3.6 m ESO telescope and thus provide much 
better spatial resolution (~ 1") than the MSX images. They 
found two compact mid-infrared sources (these are marked in 
Fig. 3) associated to G287.4779-00.5463. The brighter of these 
mid-infrared sources coincides perfectly with the disk object, 
while the fainter one can be identified with another embedded 
object with surrounding diffuse nebulosity seen in our HAWK-I 
images about 11" to the south-east of the disk object. Based on 
its J2000 coordinates of a = 10 h 44 m 32.3 s 6 = -59° 33' 14", we 
will refer to this second embedded object as YSO 32.2-14 in the 
following. 

The angular separation of the disk object and YSO 32.2-14 
is large enough that they are easily discerned in the HAWK-I, 
2MASS, Spitzer, and TIMMI2 images. For MSX and IRAS, 
however, these two objects cannot be separated and thus the 
fluxes derived for wavelengths beyond 10 /mi are contaminated 
by the contribution of the second embedded object. 



4. Morphology and spectral energy distribution of 
the Disk Object 

4. 1 . Near-Infrared Morphology 

A close-up view of the disk object as seen in our HAWK-I im- 
ages is shown in Fig. 4. The main morphological structures can 
be described as follows: 

- The general morphology is that of an approximately spheri- 
cal nebulosity with a diameter of about 4" (corresponding to 
9200 AU) and a small degree of elongation, approximately 
along the north-south direction. 

- The nebulosity is divided into a northern and a southern lobe 
(hemisphere) by a prominent dark central lane. The width 
of the dark lane is * 1.0" (2300 AU) in the /-band image, « 
0.5" (1 150 AU) in the ff-band image, and * 0.3" (690 AU) 
in the /*f s -band image. 

- While the central lane appears completely dark in the J- 
band, a faint central object becomes visible in the H- and 
,7^5 -band. 

- While the nebulosity has a rather sharp eastern edge, the 
western part is much more diffuse. This suggests a wind- 
swept morphology from some agent to the east of the ob- 
ject. We note that there are two known massive stars in this 
direction that may be responsible for this effect: The 05 star 
HD 303311, 48" east and 15" north of the disk object (at 
a projected distance of 0.56 pc), and the 03 star HD 93250 
105" east and 44" south of the disk object (at a projected 
distance of 1.26 pc). 

- The south western part of the nebulosity shows some kind of 
"tail" pointing towards the south. This structure breaks the 
otherwise high degree of symmetry of the disk object. As this 
"tail" is the only part of the infrared nebulosity that is visible 
in the optical HST images, it might be related to foreground 
nebulosity. 
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4.2. Near-infrared Photometry 

Photometric calibration of the HAWK-I images was achieved 
via bootstrapping from aperture photometry measured in the 
HAWK-I images for the star 2MASS J10442989-5933170; this 
star is about 10" to the south-west of the disk object, appears 
well isolated, and has 2MASS magnitudes of J = 14.272 ± 0.04, 
H = 13.376 + 0.04, K s = 13.015 + 0.05. 

For the bright star31.3-ll near the disk object we deter- 
mined magnitudes of / a 15.3, H as 14.6, and K s as 14.7. The 
rather bluish color of this star (H - K s = -0.1) suggests a very 
low level of extinction and may imply that this object is located 
in front of the dark cloud. It also implies that its contribution 
to the mid- and far-infrared fluxes for the disk object will most 
likely be very small and can be neglected. 

For the disk object, we determined the following measure- 
ments: The peak surface brightness in the northern / southern 
lobes is 17.85/16.97 mag/sqarcsec in the /-band, 15.03/15.01 
mag/sqarcsec in the //-band, and 12.44/ 13.13 mag/sqarcsec in 
the Zf j-band. The ratio of the integrated fluxes of the northern and 
southern lobes is 0.62 : 1 in the /-band, 0.94 : 1 in the //-band, 
and 1.47 : 1 in the A^-band. The magnitudes of the central source 
in the H- and the Z^-band (measured in circular apertures with 
2.5-pixel radii from which the average of the fluxes measured in 
two equally large apertures placed just to the east and west of 
the center in the dust lane was subtracted to account for the local 
background) are found to be K s ~ 16.1, H ~ 19.6, / > 21.8; 
note that the errors of these measurements are substantial due 
to the complicated background. If we assume that the observed 
emission is actually the central young stellar object, the observed 
color H - K s as 3.5 and the assumption that the intrinsic color of 
the central star should be (H-K s )o < 0. 1 lead to an extinction es- 
timate of Ay > 54 mag. However, we note that the central patch 
in the images may also just be reflected light from the innermost 
parts of the circumstellar disk; in that case, the true extinction 
and optical depth to the central source would be much higher. 

4.3. Morphology at other wavelengths 

In all optical HST images, star31.3-ll is clearly visible, but the 
disk object cannot be seen. Only the "tail" at the south-western 
edge of the disk object is visible in the F555W and F850LP im- 
age, suggesting that it may represent reflected light from a fore- 
ground cloud that is unrelated to the disk object. 

In Fig. 5 we compare the Spitzer IRAC images to our 
HAWK-I near-infrared images and the HST F555W image. 
Although the IRAC images do not have enough angular reso- 
lution to show the structure of the disk object, it can be seen 
that the optically invisible disk object strongly brightens with 
increasing wavelength. The same is true for YSO 32.2-14. 

4.4. Construction of the spectral energy distribution 

We assembled the spectral energy distribution of the disk ob- 
ject from the images and literature data discussed above. For 
our HAWK-I data we performed aperture photometry using an 
aperture diameter of about 5", excluding star 31.3-11. For the 
Spitzer IRAC data an aperture diameter of 7" was used. Due to 
the lower angular resolution of the IRAC data, the disk object 
and star 31. 3-11 cannot be separated. However, as the ZT,-band 
flux of star 31.3-11 is only 8.8 x 10~ 4 Jy, i.e. 40 times lower than 
that of the disk object and since star 31.3-11 shows quite blue 
colors, we are confident that its contribution to the emission in 
the IRAC bands can safely be neglected. The second embedded 



Table 1. Spectral energy distribution of the disk object. References: 
(1) This work; (2) Mottram et al. (2007); (3) Rathbome et al. (2002). 



A 




aperture 


Instrument 
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[pro] 


[Jy] 


diameter 






1.26 


0.0013 ± 0.0003 


5" 


HAWK-I 


1 


1.62 


0.0099 ± 0.0015 


5" 


HAWK-I 


1 


2.15 


0.0351 ±0.0045 


5" 


HAWK-I 


1 


3.6 


0.113 ±0.023 


7" 


Spitzer 


1 


4.5 


0.237 ± 0.047 


7" 


Spitzer 


1 


5.8 


0.443 ± 0.089 


7" 


Spitzer 


1 


8.0 


0.760 ± 0.228 


7" 


Spitzer 


1 


10.4 


0.35 ±0.01 


1.2" 


TIMMI2 


2 


8.28 


as 5.52 




MSX 


3 


12.13 


as 10.1 




MSX 


3 


14.65 


* 7.3 




MSX 


3 


21.3 


as 28.4 




MSX 


3 


25 


~ 50 




IRAS 


3 


60 


~ 300 




IRAS 


3 


100 


~ 233 




IRAS 


3 


870 


< 2.5 


36" 


LABOCA 


1 



source YSO 32.2-14 is clearly separated from the disk object in 
the IRAC images (see Fig. 5) and does not contaminate our pho- 
tometry. 

Although the Spitzer MIPS maps also have sufficient angular 
resolution to clearly separate the disk object from YSO 32.2-14, 
both objects are heavily saturated and therefore no quantitative 
photometric analysis is possible. In order to check the reliabil- 
ity of the mid- to far-infrared fluxes determined by Rathborne et 
al. (2002) from the MSX data, we analyzed and compared the 
Spitzer MIPS 24 pm image and the MSX E-band (21.3 pm) im- 
age. The MIPS image shows a bright, roughly circular nebulos- 
ity with a diameter of about 50" with three embedded compact, 
point-like sources. The brightest of these point-like sources is 
the disk object, the second brightest is YSO 32.2-14, and the third 
one is the 2MASS source 10443341-5933189 (the bright object 
near the lower left corner of Fig. 3). In order to determine the 
24 /mi flux of the nebula, we used an aperture with diameter 54" 
and excluded circular regions with diameters of 11" centered 
on each of the three point-like sources. Subtracting the large- 
scale background level estimated from several nearby regions, 
we found a flux of as 46 Jy for the nebula. In the MSX E-band 
image, the cloud is also clearly visible, but the point-like sources 
are not resolved. The total flux, measured again in a 54" diam- 
eter aperture, is as 95 Jy; this is the sum of the emission from 
the cloud and the embedded point-like sources. Subtracting the 
as 46 Jy from the nebula as estimated above from the Spitzer 
map, leaves a total flux of as 49 Jy for the sum of the three 
point-like sources. Since the disk object is clearly brighter in the 
Spitzer maps than the other two point-like sources, the flux of 
28.4 Jy determined by Rathborne et al. (2002) for the dominant 
point-like source (i.e. the disk object) is a very reasonable value. 

We also performed aperture photometry in our LABOCA 
sub-mm map. In order to approximately isolate the contribution 
of the disk object from the surrounding cloud, we measured the 
flux in an aperture with a radius of 18" (corresponding to the 
FWHM of the beam) and subtracted the scaled background flux 
estimated from the median flux level in an annulus with 1 . 1 and 
1.5 times the aperture radius. This yielded a flux of 2.5 Jy. As 
the rather large aperture partly includes YSO 32.2-14, this value 
should be regarded as an upper limit to the sub-mm flux from 
the disk object. 
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Fig. 5. The disk object and its surroundings seen at different wavelengths. The upper left panel shows the optical HST image taken through the 
filter F555; the position of the (optically invisible) disk object is marked by the white circle. The next three panels show the HAWK-I images, and 
the lower row shows the Spitzer IRAC images. 



Table 2. Spectral energy distribution of the source YSO 32.2-14 



A 


F v 


aperture 


Instrument 


[fim] 


[Jy] 


diameter 




1.26 


0.00053 +0.00011 


2" 


HAWK-I 


1.62 


0.0024 ± 0.0036 


2" 


HAWK-I 


2.15 


0.0089 +0.0013 


2" 


HAWK-I 


3.6 


0.053 ±0.011 


5" 


Spitzer 


4.5 


0.105 ± 0.021 


5" 


Spitzer 


5.8 


0.219 ± 0.044 


5" 


Spitzer 


8.0 


0.390 ±0.117 


5" 


Spitzer 


10.4 


0.24 ± 0.01 


1.2" 


TIMMI2 



The SED of the disk object (see Fig. 6) rises steeply from 
1.26 /mi to 8 /im. The spectral index measured between 2.15 /mi 
and 8 /mi is a — 0.878; the object thus qualifies as a "class 
I source", i.e. a very young object that is embedded in large 
amounts of circumstellar material. 

We note again that the reliability of the fluxes is quite dif- 
ferent in the different wavelength ranges. The near- and mid- 
infrared fluxes derived from the HAWK-I and Spitzer data as 
well as the TIMMI2 flux are highly reliable and not contami- 
nated by other infrared sources. The MSX and IRAS fluxes for 
the object derived by Rathborne et al. (2002) are less reliable due 
to the uncertain local background correction and may be contam- 
inated by emission from YSO 32.2-14. The LABOCA flux is also 
uncertain since it may contain significant amounts of large-scale 
nebula emission and is therefore regarded as an upper limit in 
our modeling described below. 

A final caveat comes from the fact that the individual points 
of the SED were not obtained simultaneously. It is well es- 
tablished that many young stellar objects show considerable 
infrared variability by factors of up to ~ 2 or even more, 



on timescales between days and decades; these variations are 
thought to be caused by changes in the structure of the inner 
disks (e.g., Flaherty & Muzerolle 2010). This implies that flux 
variations by about a factor of 2 could be present and change the 
shape of the observed SED. 



5. Radiative transfer modeling 

5. 1 . Modeling of the spectral energy distribution with the 
Robitaille models 

The pre-computed set of SED models provided by Robitaille et 
al. (2007) are now used as a standard tool to estimate stellar and 
circumstellar parameters of YSOs from the observed SEDs. We 
therefore choose to fit the observed spectral energy distribution 
of the disk source with these models as a first step of our anal- 
ysis. We note that this kind of analysis ignores the morphologi- 
cal information we have available from the near-infrared images. 
However, since most current studies of YSO circumstellar struc- 
tures are based on such SED fits without spatially resolved image 
information, this modeling can tell us what conclusions would 
be drawn if the direct spatial morphology information from our 
HAWK-I images were not available. 

Using the SED fluxes as listed in Tab. 1 and fixing the dis- 
tance to 2.3 kpc, the best fit (x 2 lv = 1.89) is found for a model 
of a young star surrounded by an envelope with the following 
parameters: M, = 13.3 M Q> L„ = 4190 L Q , M env = 42 M Q , 
Rem = 5.3 x 10 4 AU, inclination = 69.5°, and foreground extinc- 
tion of Ay = 5.59 mag. 

In the 3 1 further models with x 2 / v ~ X 2 / v (best) < 2, the 
parameters are in the following ranges: M, = 8.11... 14.3 M , 
U = 1320... 4950 L , M env = 4.8... 363 M Q , inclina- 
tion = 18.2° ... 87.1°, and foreground extinction of Ay = 
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Fig. 6. The 32 best fitting Robitaille SED models compared to the mea- 
sured fluxes of the disk object. The best fit is shown as a solid line; the 
dashed line shows the stellar spectrum of the best fit model with the 
effect of the foreground extinction of A v = 5.59 mag. The grey lines 
show the other good models as discussed in the text. 



1.24. . .20.9 mag. Nineteen of these models have (in addition 
to the envelopes) disks with masses M disk = 8.7 x 1(T 3 ...4.3 x 
lfT 1 M Q . 

We also performed SED fits excluding the near-infrared 
fluxes (that are contaminated by scattered light). We found that 
this produces only very marginal changes of the fitting results. 
The best model is the same as the second best model in the 
full SED fit. The parameter ranges for the best 3 1 models are: 
M* = 8.28... 15.6 M@, U = 1550 . . .65OOL , M enve iope = 
48.2... 251 M , inclination = 18.2° ... 87.1°, and foreground 
extinction of A v — ... 35 mag. 

The rather wide ranges for the stellar and circumstellar pa- 
rameters demonstrate the large ambiguities of the fits based 
on the SED only. If we make use of the information deduced 
from the observed morphology and restrict the selection to 
models with a nearly edge-on viewing angle (i > 75°), we 
find six different models with the following parameter ranges: 
M« = 8.45... 13.2 M e , L, = 2140 .. .4750 L Q , M envelope = 
26...93M , ^envelope = 4.4 x 10 4 . . . 1 x 10 5 AU. Two of 
these models also have a circumstellar disk; one with a mass 
of Mdisk = 0.027 M and radius = 122 AU, and the other 
one with M disk = 0.033 M e and radius 7? disk = 8 AU. 

5.2. 2D radiative transfer modeling of the SED and 
morphology 

In the second step of our analysis, we performed a more detailed 
radiative transfer modeling of the object in which both, the SED 
and the near-infrared morphology, are properly taken into ac- 
count. The direct spatial information from the NIR images can 
resolve many of the ambiguities of the models that are based on 
SED fits only, and provide a much more physically meaningful 
model. 

We used the 2D radiation transfer code RADMC 
(Dullemond & Dominik 2004), which solves the continuum 
radiative transfer problem in a passive irradiated axisymmetric 
dusty disk and envelope around a central illuminating star by 
means of a Monte Carlo algorithm. After the self-consistent 
computation of the dust temperatures throughout the disk, the 



SED and images at different wavelengths are computed with a 
ray-tracing procedure. The density model for the circumstellar 
disk and envelope we have used in our modeling is defined in 
cylindrical coordinates (r, z) by the following formulae: 
The density of the disk is given by 



Pdisk - Pdisk.O 





z 2 - 


(k) exp 


~2h? 



(1) 



where ro is set to the stellar radius and h is the vertical pres- 
sure scale height. The parameter h can increase with the radial 
distance from the center, resulting in a flaring of the disk: 



(2) 



h(r) = ho - 
Vo 



This disk is embedded in a spherical envelope with power- 
law density gradients: 



Penv — Penv.O 



and 



Penv — Penv.O 



Re 



Re 



for r < R e 



for r > R P 



(3) 



(4) 



For our computations, we used a grid size of R = 5000 AU, 
130 radial grid points, 50 angular grid points (with 10 extra an- 
gular points to refine the grid near the equator), a gas-to-dust 
mass ratio of 100:1, and 10 5 photon packets for the Monte Carlo 
simulations of each model. For the dust opacities we used the 
updated model of Draine & Lee (1984) 1 with 75% silicate and 
25% graphite. We assumed a canonical grain size distribution 
n(a) oc a -3 - 5 with a minimum grain size of 0.005 pm. For the 
maximum grain size we used a density dependent value ranging 
from a max = 0.25 pm for p < 10~ 17 g/cm 2 up to a max = 10 pm 
for p > 10~ 13 g/cm 2 . 

Our density model is described by eleven free parameters: 
(1) Stellar luminosity L*; (2) Stellar temperature T„\ (3) Disk 
density power law a; (4) Disk flaring parameter /3; (5) Disk ver- 
tical scale height ho; (6) Disk density paisk.o; (7) Inner envelope 
density power law exponent y; (8) Outer envelope density power 
law exponent 6; (9) Envelope characteristic radius R &m ; (10) 
Envelope density p en v,o; (H) Inclination ;. Furthermore, as the 
object is obviously located inside a dark cloud, we also consid- 
ered the effects of foreground extinction. 

A general problem of such a radiative transfer modeling is 
that the high dimensionality and the complicated topology of the 
parameter space make a search for the "best model" very diffi- 
cult in practice. Simple scanning of the parameter space is not 
feasible: for example, even an extremely coarse discretization of 
only 5 different values for each of the 1 1 parameters would al- 
ready require the computation (and evaluation) of about 50 mil- 
lion different models. 

Another problem is how to evaluate the fit-quality of a spe- 
cific model. For the comparison of the model SED to the ob- 
served SED, a simple x 1 analysis is easy to implement. A 
quantitative evaluation of the model images, however, is not so 
straightforward. First attempts to compare the individual pixel 
values in the model images to those in the observed images were 
not successful. Instead, we performed a quantitative assessment 
of the most important morphological features in the images. This 



1 see http://www.astro.princeton.edu/ draine/dust/ 
dust.diel.html 
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Fig. 7. SED (above) and near-infrared model images (below) computed 
with RADMC. The solid dots in the SED plot show the HAWK-I and 
Spitzer data points, the open squares the MSX data, the thick X symbol 
the TIMMI2 measurement, the open diamonds the IRAS fluxes, and the 
open square the LABOCA flux. A square-root intensity scaling is used 
for the model images. 



was implemented by computing 2D discrete cosine transforma- 
tions (DCTs) of the model images and comparing 192 DCT co- 
efficients per model image to those derived from the observed 
images. In this way we made sure that model images classified 
as "good fits" showed a central dark lane and a roundish nebu- 
losity above and below this lane. 

In the first part of our modeling, we computed several thou- 
sand models to explore the main effects of the individual pa- 
rameters on the fit quality. Based on these initial results, we 
then implemented a "genetic algorithm", in which small random 
changes of some of the parameters of good models are intro- 
duced in order to find an even better model. In total, more than 
150000 models were computed. 

It was very easy to find models that reproduce the SED very 
well. However, the images of these models deviate consider- 
ably from the observed images: they generally produced a much 
thicker dark lane, and often the shape of the upper and lower re- 
flection nebulosity was more cone-like than hemisphere-like (as 
observed). On the other hand, we also found several models that 
reproduced the near-infrared morphology quite well; these mod- 
els, however, could not reproduce the shape of the observed SED 
between 2 /mi and ~ 20 /mi. 

We finally found several models that reproduce both, SED 
and images, in an acceptable way (although none of these 
models could simultaneously reproduce the SED and the near- 
infrared morphology perfectly). One of these "reasonably good" 
models is shown in Fig. 7. The model parameters are: 

- Stellar Luminosity: L, = 12 000L o 

- Stellar temperature: T t = 22 000 K 

- Disk density power law: a = 2.25 



- Disk flaring parameter: J3 = 1 .25 

- Disk vertical scale height: ho/ro = 0.2 

- Total disk mass: M d i s k = 1.78 M Q 

- Inner envelope density power law exponent: y — 1 .3 

- Outer envelope density power law exponent: 5 = -4.5 

- Envelope characteristic radius: R env = 1500 AU 

- Total envelope mass: M env = 0.16 M 

- Inclination: i — 87° 

The model assumes a foreground (i.e. interstellar cloud) extinc- 
tion of Ay = 2.5 mag. Comparison of the stellar luminosity and 
temperature to the pre-main sequence models for massive stars 
from Bernasconi & Maeder (1996) suggests a stellar mass of 
« 12 M Q . Such an object should evolve towards a ZAMS star of 
spectral type w B 1 . 

Similarly good models can be found for the following range 
of the main parameters: L t a 10000 ... 15 000 L , M d i s k ~ 
1.8 . . .4.7 M , and M env « 0.16 . . .0.5 M G . The range of val- 
ues for the vertical disk scale height in the good models is 
ho/ro — 0.1 . . .0.2; this agrees well with the typical scale heights 
found for circumstellar disks around young stars (see Dullemond 
& Dominik 2004) and clearly shows that the spatial configu- 
ration of the circumstellar material is that of a geometrically 
thin disk, and not a thick torus (that would be characterized by 
ho/ro > 0.5). The luminosity range suggests that the most likely 
mass of the central object is in the range « 10 ... 15 M Q . 

One fundamental difference between this model (for SED 
and near-IR image morphology) and the Robitaille models (for 
the SED only) is the presence of a massive and optically thick 
circumstellar disk. Such a disk is not required for a good fit to 
the SED alone, but its presence is clearly established from the 
near-infrared images. Another notable difference is seen in the 
values of the stellar luminosity, which is considerably higher 
in our RADMC models compared to the Robitaille SED mod- 
els. This difference is related to the so-called "flashlight effect" 
first discussed by Yorke & Bodenheimer (1999). It results from 
the highly non-isotropic distribution of radiative flux that oc- 
curs when a dense, optically thick and flat circumstellar disk is 
present around an illuminating central source. At shorter wave- 
lengths (for which the disk is optically thick), the radiation field 
is strongly concentrated toward the polar directions. In such a 
situation, the total observed outgoing radiation flux is a strong 
function of the viewing angle; the integrated observed flux for 
nearly edge-on inclinations ; « 90° can be considerably smaller 
than the integrated observed flux measured from a more pole- 
on view position (see Yorke & Sonnhalter 2002). In the case of 
the Carina disk object, the presence of the dark central lane di- 
rectly shows that much of the stellar light is removed from our 
viewing direction and scattered into directions perpendicular to 
the line-of-sight. Therefore, a high stellar luminosity is required 
to reproduce the observed flux seen in our edge-on viewing di- 
rection. The Robitaille models (that use only the SED as input), 
on the other hand, have no such optically thick disk around the 
star and thus assume a much more isotropic radiation field. As a 
consequence, these models underestimate the luminosity of the 
object. 

It has to be emphasized that the detection and the 2D char- 
acterization of the disk object was only possible due to the su- 
perb image quality of the HAWK-I data. In order to check how 
strongly our modeling results depend on the angular resolution, 
we applied Gaussian filters with different FWHM to the images 
to see how the observed morphology changes. Degrading the in- 
trinsic resolution of the images 0.4" in the i^-band) to a more 
typical value of 0.7" is already sufficient to reduce the contrast 



Preibisch et al.: A large massive circumstellar disk surrounding a high-mass young stellar object in the Carina Nebula 



9 



between the dark lane and the bright lobes to less than 10%. This 
implies that under "average" seeing conditions the presence of 
the dark lane (and thus the nature of the object) would no longer 
be clear. Since the dark lane is the strongest observational con- 
straint for our 2D modeling, the existence of a disk would no 
longer be established. 



6. Summary and Conclusions 

Our near-infrared survey of the Carina Nebula led to the 
serendipitous discovery of an interesting embedded young stel- 
lar object that is surrounded by a circumstellar disk and a very 
large envelope. The presence of the prominent envelope suggests 
a very early evolutionary stage of the young stellar object. The 
most notable characteristics of this object derived from our mod- 
eling can be summarized as follows: 

The central young stellar object has a high estimated mass of 
about 12 M ; this object is thus one of the very few known cases 
of massive young stellar objects that are clearly surrounded by a 
circumstellar disk. 

Another outstanding feature is the very large physical size of 
the disk, which has a radius of 2700 AU. This is thus much larger 
than most disks found around lower-mass young stellar objects 
(that have typical radii of a at most a few hundred AU). 

Finally, our modeling suggests a quite high mass of about 
2 M for the circumstellar disk, and an additional envelope mass 
of about 0.2 M . The total mass of the circumstellar material 
(within $ 10 000 AU) seems to be of the order ~ 2 - 3 M Q . This 
is a remarkably high value for the circumstellar material around 
an individual young stellar object. Also, the ratio between cir- 
cumstellar mass and stellar mass is high, ~ 20%. This is much 
higher than typically found for lower-mass young stellar objects. 

Based on these properties, the Carina disk object clearly rep- 
resents one of the most interesting circumstellar disk detections 
made so far. Further observations are planned to study this ob- 
ject in more detail. An important point is the uncertainty of 
the current flux determinations at wavelengths above ~ 10 //m, 
where the angular resolution of the existing data sets does not 
allow an accurate separation of the disk object from neighboring 
embedded objects and the surrounding diffuse cloud emission. 
Although we have tried to account for these effects, it is still pos- 
sible that the SED fluxes are contaminated by the surrounding 
nebulosity (as in most other far-infrared observations of embed- 
ded objects). Therefore we cannot exclude the possibility that 
the fluxes from the object are actually somewhat lower, which 
could lead to smaller values for the disk mass and the mass of 
the YSO. 

The radiative transfer modeling demonstrates the essential 
value of direct image information. As mentioned above, the de- 
tection and characterization of the disk was only possible due 
to the very good angular resolution of our HAWK-I images and 
would not have been feasible under more typical seeing condi- 
tions. This emphasizes (once again) the importance of angular 
resolution in the characterization of circumstellar matter around 
young stellar objects. Images with even higher angular resolu- 
tion (e.g., by adaptive optics observations) would provide very 
interesting new information about the width and the shape of 
the disk shadow. In particular, the chromaticity of the brightness 
distribution along cuts perpendicular to the disk plane could pro- 
vide crucial new constraints for a more detailed radiative transfer 
modeling, and yield new insight into the physical parameters and 
dust properties of the circumstellar material around this massive 
protostar. 
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